	set more off

use Data/Estimation_sample.dta, clear
	
	gen nexis_buildings  = ROOFME/PROOFME
		
		
	gen roofratio = ROOFME/ROOFTI
	gen roofratiobuildings = roofratio*nexis_buildings
		
	gen roofrelativeIV = (roofratio*nexis_buildings/[moduleprice])/1000

	

	
	rlasso green_p cuminst i.quarter median_income_interp median_mortgage_interp median_rent_interp PBachelor_interp PFullTime_interp median_age_interp median_rooms_interp PSeparate_interp  PRented_interp [aweight =NCustomer] , cluster(postcode) fe pnotpen(cuminst )
	
	eststo lasso1
	
	
	rlasso green_p roofrelative i.quarter median_income_interp median_mortgage_interp median_rent_interp PBachelor_interp PFullTime_interp median_age_interp median_rooms_interp PSeparate_interp  PRented_interp [aweight =NCustomer] , cluster(postcode) fe pnotpen(roofrelative )
	
	eststo lasso2
